Multi-Partite Entanglement Inequalities via Spin Vector Geometry 
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We introduce inequalities for multi-partite entanglement, derived from the geometry of spin vectors. The 
criteria are constructed iteratively from cross and dot products between the spins of individual subsystems, each 
of which may have arbitrary dimension. For qubit ensembles the maximum violation for our inequalities is 
larger than that for the Mermin-Klyshko Bell inequalities, and the maximally violating states are different from 
Greenberger-Horne-Zeilinger states. Our inequalities are violated by certain bound entangled states for which 
no Bell-type violation has yet been found. 
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Entanglement is one of the most mysterious features of 
quantum physics and a key ingredient in the science of quan- 
tum information. While initial research was focussed on bi- 
partite entanglement, multipartite entanglement has attracted 
increasing attention, because it was realized that multipartite 
entangled states can exhibit qualitatively different features [ 1 ] . 
Multi-partite entangled states are also important for most ap- 
plications envisaged in quantum information, such as quan- 
tum computation [2]. Recently, multi-partite entanglement 
has been studied using multi-partite Bell inequalities [3, 4], 
the partial transposition criterion [6], and a variety of other 
methods [7, 8]. 

Here we develop a new approach based on the geometry of 
spin vectors. Our results are relevant for any system where 
operators analogous to spin can be defined and measured, e.g. 
for multi-mode light fields [9] and Bose-Einstein condensates 
[10]. In any finite-dimensional Hilbert space one can always 
define operators { J x , J y , J z } = J that satisfy the commuta- 
tion relations of angular momentum, [J u , J v ] = ie uvw J w for 
u, v,w £ {x,y,z}. The transformations U(a) = e ta ' J gen- 
erated by these operators, where a is a numerical vector, form 
a representation of the group SU(2). Applying U{a) to the 
quantum states leads to an SO(3) rotation R(a) of the vec- 
tor (J) = {{J x ), {Jy), {Jz)} of spin expectation values. The 
vector ol gives both axis and angle of rotation. 

Dot Product: Our inequalities involve the expectation val- 
ues of operators that are constructed from the spins of sub- 
systems via the cross and dot product. We will first illustrate 
the principle with the simplest example, the dot product be- 
tween two spins of magnitude j\ and j 2 . For a product state 
\1p12) = \<Px) \Xz) the expectation value of the dot product is 



changing its norm. One has \{Jz }| < jx, because jx and 
—jx are the eigenvalues of Ji with the largest modulus. As 
a consequence, we have for every product state 

\{DW)\ = \{J^).jW)\/(j x j 2 )<l. 



(2) 



The bound of Eq. (2) can be extended straightforwardly to 
separable states given the following triangle inequality (the 
expectation value for product states has suffix i/j): 



(3) 



Note that the choice of co-ordinate system for each spin is 
arbitrary. 

What is the maximal violation of Eq. (2) for entangled 
states? Noting the relation (jW.j( 2 )) = ±(J 2 - Jt 1 ) 2 - 
j( 2 ) 2 ), with J = J^ 1 ) + j( 2 \ indicates that the eigenstates 
of (J^.J^ 2 )) are also eigenstates of J 2 , corresponding to 
j(j + 1) where j = j 1 + j 2 - A, and A e JN < 2Min[j" 1 , j 2 ] . 
Assuming ji < j 2 and that subspaces (1) and (2) have fixed 
dimension gives: 



1 + ^ — ({J^.J 1 ^)) <1 (4) 
32 J 3x32 V / 



From this result it is seen that the maximum absolute value is 
3, provided by the 2-qubit singlet, j\ = j 2 = 1 /2. 

Cross Product and Multipartite Inequalities: In full anal- 
ogy with the dot product, the expectation value of the cross 
product in a product state is the cross product of the spin ex- 
pectation value vectors for individual systems, 



{i>X2\J^ 2) X jW|^ 12 ) = (X2|J I2) |X2> X {falfV 



(5) 



(feljW • J(%12) = W>l|J (1) |0l) • (X2|J (2) |X2), CD 

which is the scalar product of two vectors {J 1 - 1 ') = 
{4>i\J {1) \(t>i) and (J( 2 )) = (X2|J (2) |X2)- The modulus of 
the scalar product is bounded by the product of the norms of 
the two vectors. Furthermore the norm ||(J^)|| cannot ex- 
ceed j\. This can be seen by noting that by a rotation the 
vector ( J^>) can always be brought to a form where only one 
of its components, say (ji ), is different from zero, without and again the generalization to separable states is immediate. 



As the norm of the cross product of two vectors is bounded 
from above by the product of the their norms, we find for 
product states 



\\{CW)\\ = \\{Wx^)\\/(j 1 j 2 ) <1, 



(6) 
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(2,1/2) 
(3,1/2) 
(4, 1/2) 
(5,1/2) 
(6,1/2) 


2 

2y/2 (w 2.828) 
2^ (« 4.899) 
2vTi (« 7.483) 
w 12.144 


3 

2^ (« 3.464) 
4^ (« 6.928) 
4^6 (« 9.798) 
w 16.971 


(2,1) 
(3,1) 
(4,1) 
(5,1) 
(6,1) 


V2 (« 1.414) 
v/3 (w 1.7321) 
\/ 3 + a/5 (w 2.288) 
« 2.840 
« 3.731 


2 

v/3 (« 1.7321) 
2^ (» 2.828) 
3 

w 4.472 



TABLE I: Some maximal violations for the inequalities of Eq. (7) 
and Eq. (9), for entangled states of N qubits (J = 1/2) and qutrits 
(J = 1). The largest violation we were able to find numerically was 
Max(L» (11) ) « 152.691, for 11 qubits. D (11) has 2 11 eigenvalues 
with some degeneracy, see Eq. (11). 

By iterating the cross product operation it is possible to de- 
rive bounds for multipartite systems. For a fully separable 
state of N spins ji,j2,ja,—,jN one derives the following 
bound, in analogy with Eq. (6): 

||(CW)||<1 (7) 

where we use the notation: 

CW = JW x (J^-D x ...( x jV))Khh.~3N). (8) 

Considering the dot product between a single spin and a vector 
constructed like in Eq. (8) one can also derive a bound 

\(D^)\<1. (9) 

for fully separable states, where 

#» = ^.(^-^..x^xjW))/^..^). (10) 

We have investigated how strongly the bounds of Eq. (7) and 
Eq. (9) can be violated by entangled states, TABLE I. For 
J)(") this involves finding its largest eigenvalue. Maximum 
values of || 1 1 may be found by studying the greatest ex- 
pectation value for any of the vector components of C^ N \ 
This is because joint identical SU(2) transformations on all 
spins (i.e. each with the same a) correspond to simple rota- 
tions of the vector (C^ N '). It can thus always be brought to a 
standard form where only one of its components e.g. Cz 
is non-zero, without changing its norm. The maximum of 
the norm is therefore the largest eigenvalue of Cz ■ Upper 
bounds also exist for all entanglement partitions and one par- 
tition may have a range of upper bounds depending on the or- 
dering of sub-systems (1), (2), . . . , (AT) in the directed prod- 
ucts flW and & N \ see TABLE II. 

Symmetry and Eigenstates: In addition to the numerical re- 
sults, some insight into the structure of the eigenstates may 
be gained from symmetry considerations. (£) W) and (C^) 



transform like a scalar and a vector respectively under all iden- 
tical joint rotations R(a)® N , or equivalently under all identi- 
cal local SU(2) transformations of the state. Furthermore, 
D( N > is anti-symmetric under the permutation (1) *-> (2) be- 
cause (jW x j( 2 )) = -(J( 2 > x J«). Operator DW must 
therefore have the following highly symmetric structure: 

£w = £m4 12) -i4 21 >) (id 

v 

i.e. a weighted sum of projectors Hp onto spaces associ- 
ated with irreducible matrix representations T> of SU (2). All 
the projectors are orthogonal, IIx>IIxy = Sx>v^v and they 
each project into a (2 J + 1) dimensional space of overall spin 
</(£>)■ States in each representation V have a shared permu- 
tation symmetry of all the A~ particles; e.g. the highest spin 
representation T>* , for which J = j'W + . . . + jW, is in- 
clusive of all states symmetric under all particle permutations. 
The projector IF 21 ) is formed from II^ 12 ) by exchanging par- 
ticles labelled '(1)' and '(2)', mapping either ITp 2 ^ to itself, 
in which case it vanishes from Eq. (11), or to an orthogonal 
projector of the same spin. Thus, eigenvalues of D^ N > appear 
in pairs of opposite sign ±(/i-p — Mx>')> w ith the multiplicity 
of pairs 2J(T>) + 1, To give a concrete example, eigenstate 
|^4) maximally violates |(Z?( 4 ))| < 1: 

2V6|^ 4 ) = (i + V3)(| TTU) + |UTT)) 
+(1 - V3)(l ITIT) + I tlTD) - 2(| ITU) + I TUT)) (12) 

This state is one of two orthogonal J = states of 4 qubits, 
and gives L> (4) |i/> 4 ) = 4^|V> 4 )- 

Since C?W is a vector operator i.e. [C^ N \ J v ] = 
ie uvw Cw for u,v,w £ {x, y, z}, it is a spin-1 object, of 
which C z is the m = component. As a consequence and 
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Max||<C< 4 >>|| 


Max|(/J (4) )| 


Max(F (4) ) 


[1 2 3 4] 


2^6 (« 4.899) 


4^3 (w 6.928) 


2y/2 (~ 2.828) 


[1 2 3|4] 


2^2 (« 2.828) 


2^2 
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[1 2|3 4] 
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[1|2|3 4] 
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[1|2|3|4] 
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TABLE II: Magnitudes of (C (4) ) and (L> (4) ) have distinct upper 
bounds (found numerically) for entanglement partitions of 4 qubits, 
labelled '(1)' to '(4)'. No entanglement exists between qubits sep- 
arated by a vertical bar. There are 4! permutations for the directed 
products (7' 4 ' and JJ>' 4 \ resulting in a range of upper bounds for 
some partitions, the largest of which is shown above. For example, 
[1|2|3 4] has three distinct bounds for Max||(C (4) )||, namely 1, y/2 
and 2. Note that the Mermin-Klyshko operator gives identical 
upper bounds [5] for [1|2|34] and [12j 34] ; its expectation value is 
also degenerate under all particle re-orderings. The new inequali- 
ties produce a larger maximal violation. Because the new partition 
bounds are the result of a (global) numerical optimization, greater 
computing power and a more refined search may allow some of them 
to be improved. 
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in contrast to D^ N \ operator Cz is not diagonal in the spin 
basis. The Wigner-Eckart theorem [11] may be invoked to 
reveal that 

(J',m'\C^\J,m) = (J,l;m,0|J,l; J',m!) T JtJ , (13) 

where (Ji,j2', Wi,Tn2|ji,j2; j, m) are the Clebsch-Gordan 
coefficients and Tj t ji is a transition matrix element dependent 
only on J and J' . The spin-1 selection rules are the familiar 
ones of dipole radiation; (J' — J) and (ml — m) E {— 1, 0, 1}, 
with the J = J' = transition forbidden, i.e. T = 0. From 
this perspective, the scalar operator D^ N > is a spin-0 object, 
only able to couple states for J = J' and m = m' . 

The pair anti-symmetry of is also true of Cz (eigen- 
values appear in pairs ±x). An example of the states maxi- 
mally violating Eq. (7) is \4>4), 



W3|0 4 ) = 3(| TITT) - 
+\/6(| TUT) + I ITU) - I TITI) - 
+1 HIT) - I TTTI) + I IITI) - I TUT) + I Till) 



ITU)) 
ITIT)) 
I ITTT) 

(14) 



which gives Cl 4) |0 4 ) = 2V6|0 4 >, cf. TABLE II. This is one 
of four orthogonal states that give the same maximum. 

Comparison with Mermin-Klyshko Inequalities: There are 
2 independent Bell inequalities for N qubit ensembles hav- 
ing two two-valued observables per qubit or 'site', [12]. This 
set is the simplest and best understood of multi-partite Bell 
inequalities, although others may be formulated with e.g. 3 
observables per site, [13]. The 2 2 inequalities are satisfied 
by all local hidden variable theories, and all are maximally 
violated by GHZ states: 



I GHZ 



N I 



(IT) 



\<g>N 



(15) 



Note that the states which maximally violate Eq. (7) and Eq. 
(9), e.g. |04) and \ip4), are generally inequivalent to IGHZ4) 
under local unitaries. This can be proved by determining the 
Schmidt coefficients of the states for a bipartite {2, 2} cut. 

For the Bell inequalities of the type mentioned, the 
Mermin-Klyshko (MK) inequality [3] has the largest possi- 
ble violation [12], within the context of quantum mechanics. 
The MK inequality depends on operator functions F W that 
may be defined recursively: 




FIG. 1: Ratio of successive maximal violations for qubits and (inset) 
qutrits. R(N) i-> Max||((7 (iv, )ll/ Max IK < 5 (JV_1) )ll (unbroken line), 
and^ Max|(D (JV) )|/Max|(D (Ar-1) }| (dashed line). D (1) is unity, 
and = J/j. The ratio for Mermin-Klyshko inequalities is \/2 
(grey horizontal line). 



the ratio of maximal violation for N qubits to (N — 1) qubits 
is always \[2 for the MK inequalities, whereas for the inequal- 
ities of Eq. (7) the ratio is always > \/~2 but displays an 
intriguing 'see-saw' character; FIG.l gives ratios for qubits 
and qutrits, and for both Dot and Cross inequalities. That a 
larger violation is possible for the new inequalities may be at- 
tributable at least in part to having three (orthogonal) observ- 
ables per qubit, compared to two for the MK inequalities. An- 
other possible reason is that our inequalities do not exclude the 
possibility of a local hidden variable model. They are strictly 
criteria for non-separability. 

Operators C^ N > and are also unlike in that they 
are not symmetric under particle exchange; C^ N > is a 'di- 
rected' product. TABLE II shows how different particle or- 
derings results in different bounds on and . 

In terms of the numerical search needed to find a violation, 
three parameters will define 3 orthogonal directions in R 3 , 
whereas 4 variables are needed for both measurement direc- 
tions per subspace in a Bell inequality. Thus, fewer parame- 
ters are required in the optimisation of Eq. (8) and Eq. (10). 

Entanglement detection and Robustness: Consider ~ 
j(3) .(j( 2 ) x JW) j 'jxjzjz for 3 qubits. Its largest amplitude 
eigenvalue is 2^3, an associated eigenstate is 

\w 3 ) = ( I TTI) + e ia \ TIT) + e t/3 | ITT))/V3, (17) 

a state inequivalent to |GHZ 3 ) under rotations of the local 
coordinate systems. For GHZ states, numerical optimisa- 



2F (N) = F^-V^iA^+A^^+F^-^^iA^-A^) tion over all local coordinate systems gives |(D (3) )| < §^3, 



(16) 

where A^ N > and A^ are observables of the Ath qubit with 
eigenvalues ±1, and F^ = A^\ Operator i^W is obtained 
under exchange A <-> A for all the nested observables. The 
MK inequality satisfied by local hidden variable theories is 
(FW) < 1. 

For any quantum state, F^ N > has an upper bound, (F^ N ') < 
2< Ar - 1 )/ 2 i resulting in a smaller possible violation of the MK 
inequality above than that attainable for Eq. (7) and Eq. (9) in 
qubit ensembles N h-> {2, 3, 4, ... , 11}, cf. TABLE I. Also, 



i.e. the maximum possible violation is smaller than that for 
the W state by a factor of 3/4. For completely separable 
states, |(Z?( 3 ))| corresponds to the volume of a parallelop- 
iped with sides of unit length. The detection of VF-type en- 
tanglement is robust against noise: Mixed with white noise, 
(1 — ^)|W / 3)(W / 3| + v Is/8, the fraction of noise v can be as 
high as 71% and |(-D (3) )| < 1 is still violated. Substituting 
IGHZ3), violation occurs for v < 61%. The 3 qubit MK in- 
equality will only detect IGHZ3) mixed with less than 50% 
Is/8, even though it is maximally violated by such states. 
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Violation Ratio: We now show that for a given entangled 
state the maximum possible violation for the Cross criterion, 
when optimizing the choice of local coordinate systems, can- 
not be larger than the maximum violation for the Dot cri- 
terion. Consider the correlation i <G {x,y,z} defined as 
i = ( 4 N) ® C\ N ~ 1] ) /jW. One may write: 



(N) 



c 



i£{x,y,z} 



" 1) )/j (JV) =x + y + z (18) 



A correlation set {x, y, z} may be mapped into 
{x, — y, — z}, {— x, y, — z} and {— x, — y, z} by local 
unitaries (tt rotations of the TVth qubit about x,y,z axes 
respectively). Therefore Max|(-DW)| is associated with 
correlations (J 4 ® c\ N ^) all having the same sign, i.e. 
Max|(£)W)| = Max(|x| + |y| + \z\). For the Cross criterion, 
taking Max 1 1 | = Max| (C { Z N) ) |, one has: 



Max|(C( w) )|=Max 



t(JV) 
j(N) 



Max 



(19) 



(which is Maxjx + y|) because a local rotation in the TVth 



subspace transforms { J X N \ J y N ^} i— ► { — Jy JV; , Ji"" 1 } above. 
Thus for a given state, the maximum of |(ci )| for all 
choices of local coordinate systems is at most a sum of two of 
the absolute values, |a;|, \y \ and \z\, obviously upper-bounded 

by Max(|x| + \y\ + \z\) 



Max||(C (Ar) )|| < Max\(D (N) )\ 



(20) 



It is stressed that for Eq. (20) the particles or subspaces are 
considered in the same order for both and flM. 

Bound Entanglement: We consider a mixture of GHZ and 
product state projectors [6]: 



Pn = 



N 



X — ^| GHZ 7v ) (GHZyv I + \ YS Un + (21) 



Here II„ is a projector onto product state 
I T)i| T>2 - - - 1 i)r*...| T)jv, i-e. only the nth qubit is 
in the orthogonal state. Projector n fl is obtained from n„ 
by interchanging all | with J,, In [6] it was shown that p^ is 
entangled for N > 4 (it has negative partial transposition in 
some partitions) but entanglement cannot be distilled for any 
{1, (N - 1)} partition, it is 'bound' [15]. See also [7]. The 
state violates (F ) < 1 if and only if N > 8 and violates a 
Bell inequality of three dichotomic observables per qubit [13] 
for N > 7; recently it was shown to violate a 'functional' 
Bell inequality [16] for N > 6. Bell violations for lower 
N have yet to be shown. In contrast, the N = 4 state 
violates both 'Cross' and 'Dot' inequalities: A numerical 
search over local unitaries gives Max||((7( 4 ))|| w 1.09 and 
Max|(L»( 4 ))| w 1.25 for p A . This result was unexpected, 



compared with the non- violation of MK inequalities by p^. 
After all, the only entanglement in pjy is due to the GHZ state 
- this maximally violates the MK inequality. 

Summary: A simple geometric approach allows the formu- 
lation of entanglement inequalities based on the scalar and 
vector products of two spin operators. This idea was extended 
to inequalities for multiple subsystems, each of arbitrary di- 
mensionality. The maximum violation for N qubits is greater 
than that for certain classes of Bell inequalities, including the 
Mermin-Klyshko inequalities. The maximally violating states 
are not generally GHZ states; examples were given and ele- 
ments of their structure discussed. showed a high level 
of robustness in detecting both W and GHZ entanglement. 
Maximal violations for all entanglement partitions of 4 qubits 
were found numerically - none of them is degenerate (un- 
like the MK bounds), and maximum violation is for a fully 
4-entangled state. We have also shown that the new inequali- 
ties can detect bound entanglement in states for which a Bell 
violation has not been found. The work here followed a spin 
description, but may be applied to arbitrary systems for which 
{Jx, Jyi Jz} operators can be defined and measured. 

We would like to thank D. Bouwmeester for useful com- 
ments. 
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